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Strongly Coupled Quark-Gluon Plasma: 
Equation of State near T c * 
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We test the quark mass dependence implemented in the quasiparticle dispersion relations 
of our quasiparticle model for the QCD equation of state by comparing with recently available 
lattice QCD data near T c employing almost physical quark masses. In addition, we empha- 
size the capability of our model to successfully describe lattice QCD results for imaginary 
chemical potential and to analytically continue the latter to real chemical potential. 
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I. INTRODUCTION 

The equation of state of strongly interacting matter is of paramount interest for understanding 
and describing the dynamics of relativistic heavy-ion collisions and the early universe and compact 
stellar objects as well. The grand canonical potential fl depends on a set of intensive quantities 
(temperature T, various chemical potentials accumulated in /2), and parameters like quark masses, 
flavor number Nf and color number N c . For QCD as theory of hadrons and quarks and gluons 
these parameters are fixed, but it is instructive to study their impact on Q, given the complexity 
of QCD. Of particular interest is £7 near the demarcation line (phase boundary) T c (/2) of hadrons 
vs. quarks and gluons. This region is probed in heavy-ion collisions [l| and is traversed in the 
cosmic evolution [2]. 

To parameterize the equation of state near T c we employ a quasiparticle model (QPM) [j], y, [2, 
HI & 0, @| which is adjusted to available lattice QCD data. The striking simplicity of the model 
consists in the expressions for net baryon density ub and entropy density s = ^2 i=a „ si reading 
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where d q = N c Nf, d g = iV c 2 - 1, e q = 1, e fl 



V is the volume, is to be calculated consistently with Eqs. HI |2|). 
attention onto considering one independent chemical potential /j, q . 



1 and /j, g = 0. The pressure p = —ft/V, where 
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We restrict here our 
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n qyg (T,n q ) are 
the dispersion relations for the quasiparticle excitations with self-energies n 9 . g (T, (i q ) in line with 
1-loop approximations and a procedure to include nonzero quark masses according to [9( 



n, 



9:9 



q,9 



2m qig u)g !g 



20? 
q,9 



(3) 



(0) 



where m g = 0. In the quark sector, the rest masses m q contain a "true" rest mass term m q and 
an artificial "lattice" mass term t; q T introduced for calculational purposes on the lattice, to 2 , = 
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1 """ £ 2 ^ 2 - The interaction generates the self-energy contributions Cj 2 g = G 2 (a q . g T 2 + P q , g /J- 2 ) 
parametrized at fi q — via 



with a q = i, j3 q — g^-, a g = ^(2iV c + JV/) and /3 g = ^. G 2 is an effective coupling strength 



fGl_ loop (C(T)), T>T C , 

G 2 (T) = <^ (4) 

[Gtioo P (C(r c )) + 6(i-§), t<t c 

with £(T) = \(T — T S )/T C , and approaches smoothly the perturbative region at large temperatures. 
Near T c , G 2 becomes large, and the shift parameter T s regulates the coupling. Below T c , the 
coupling changes drastically its behavior. 

The model may be considered as a resummed expression for the thermodynamic potential, 
as a formal power expansion in G 2 generates an infinite series of terms, including also a term 
proportional to the plasmon term, for instance. It goes beyond a perturbative expansion scheme 
as result of the effective coupling G 2 , which also may repair possible shortcomings of the 1-loop 
inspired parameterizations of the dispersion relations. 

In this form, the model does not contain critical point (cf. [HI) or color- flavor locking effects (cf. 

El)- . 

Despite of its simplicity, the equation of state for real and imaginary chemical potential as well 
as various susceptibilities are described very well [12|, LUa] in the QPM. Here, we describe one new 
aspect of our model, namely a naive chiral extrapolation of the equation of state. We emphasize 
also the capability to extrapolate lattice QCD results from imaginary /j, q to real /i q . We focus on 
the region T ~ T c , where interaction effects are strong. It is the region of the strongly coupled 
quark-gluon plasma [14 , lla . Il6| presently explored experimentally at RHIC and in near future at 
LHC and later on at FAIR. 



II. CHIRAL EXTRAPOLATION 

The QPM parametrization of the lattice QCD results [13| for Nj = 2 + 1 with fairly large quark 
masses, i.e. £, u ,d — 0.4 and £ s = 1, was already presented in [l| for the scaled pressure p/T 4 at 
n B = 0. The QPM parameters read A = 7.8, T s /T c = 0.8, b = 347 for m { ^ ds = 0. 

Now we try to accomplish a chiral extrapolation by means of Eq. ([3]). Neglecting naively a 
conceivable dependence of A, T s /T c and b on m q , the extrapolated results are exhibited in the left 
panel of Figure Q] (dashed curve) for £ Mj d = 0.015 and £ s = 0.15, corresponding to the set-up in 
J18j |. These lattice QCD results [lj| (squares) are astonishingly well reproduced, however, by the 
price of changing the pressure integration constant, B(T C ), which needs to be readjusted in order 
to generate the small pressure below T c . We note that putting £ u ,d.s = in the strikt chiral limit 
does not change noticeably the dashed curve on the scale exhibited in Figure 1. 

Another important test [19|] of the suitability of our hydrodynamic EoS is the interaction measure 
(e — 3p), where e denotes the energy density. As shown in the right panel of Figure[T] the QPM for 
almost physical quark masses (i.e., £, Ut d = 0.015 and £ s = 0.15) faithfully reproduces corresponding 
lattice QCD data [18| of (e — 3p)/T 4 . The peak, which is related to the softest point in the QCD 
equation of state, is located at T/T c = 1.08. For larger temperatures, the interaction measure 
approaches logarithmically zero according to the temperature dependence in the effective coupling 
G 2 , though, is close to the conformal limit e = 3p already for T/T c > 10. 
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FIG. 1: Exploration of the quark mass dependence in the quasiparticle dispersion re- 
lations. Left: Circles and squares exhibit lattice QCD results for the scaled pressure 
p{T)/T* for Nf = 2+ 1 at n B = for m u>d = 0.4 T, m s = T [H and almost physical 
quark masses yj| , respectively. The black solid curve shows the QPM parametriza- 
tion [1| of lattice QCD data [13]. The red dashed curve exhibits the QPM result 
when changing the quark masses to m u ,d = 0.015 T, m 3 — 10 m u .d, corresponding to 
[1SJ |. leaving the parameters in G unchanged but rendering B(T C )/T C from 0.54 to 
0.76. Right: Comparison of QPM with lattice QCD results for the scaled interaction 
measure (e — 3p)/T 4 (line codes and symbols as in the left panel). 



III. IMAGINARY CHEMICAL POTENTIAL 



In this section, QCD thermodynamics for Nf — 4 at nonzero imaginary chemical potential 
/iq = ifj,i is considered. The net quark number density n q — 3 ub reads in the QPM [131 ] 



n q (T,ifii) = i-f 
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which is purely imaginary and an odd function in /x^. In Figure [2] (left panel), QPM results for 
—irigjT 3 , employing the parametrization from [13| ]. are compared with the lattice QCD results 
[2fJ, |21(. In particular, the pronounced bending for T = 1.1 T c close to \x c jT = 7r/3, which signals 
the onset of the first-order Roberge- Weiss transition [23[ , is accurately reproduced representing a 
sensible test of the QPM at nonzero Ub- In the QPM, the bending is driven by the increasing /u, 
dependence in the quasiparticle dispersion relation io q close to /i c . For temperatures T > 1.5 T c 
we observe an independence of n 9 /T 3 considered as a function of Hi/T on the explicit value of T. 
This independence follows from Eq. ([5]) as long as uj q is approximately independent of fii. 

The phase boundary T c {^jlb) for Nf — 4 is shown in the right panel of Figure [2] for imaginary as 
well as for real baryo-chemical potential fis = 3/i g , which is accessible by analytic continuation. 
Within the QPM, an estimate for T c (/j,b) follows as self-consistent solution of a partial differential 
equation based on Maxwell's relation emerging at T c (fis = 0) (cf. [13|]). T c (/j,b) for imaginary 
[is and the first Roberge- Weiss transition line cross each other at T*/T = 1.112 and (up) 2 /T? = 
-12.214 being close to the lattice QCD results T* /T = 1.095 and (^) 2 /T c 2 = -11.834 UlSl]. 

The features of the phase diagram in the imaginary chemical potential sector can be described as 
follows. Dashed curves represent the analytic sections of the first two Roberge- Weiss transition lines 
H B /Tc = — T 2 7r 2 (2fc+1) 2 /T c 2 (here k = 1, 2) turning into first-order transitions (solid section) while 
the dotted curve shows the first Z$ center symmetry line /J. 2 3 /T c 2 = — 4T 2 7r 2 /T c 2 . The Roberge- 
Weiss periodicity [23| implies that sectors between fJ,i/T — 2-7rfc/3 and 27r(fc + l)/3 (sector II for 
k = 1) are repeated copies of sector I between fXi/T = and 27r/3. The subsector between 
Hi/T — 7r/3 and 2-7r/3 (sector I") is a reflected copy of the subsector between fXi/T — and 7r/3 
(sector I') mirrored at the first Roberge- Weiss transition line. 
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FIG. 2: Left: Comparison of QPM [13[ (solid curves) with continuum estimated 
lattice QCD results [201 . |21| | (symbols) for the imaginary part of the scaled net 
quark number density rtq/T 3 for Nf = 4 as a function of fXi/T for temperatures 
T — 1.1, 1.5, 2.5, 3.5 T c (diamonds, circles, squares and triangles, respectively). The 
discontinuity in n q /T 3 takes place at the Roberge- Weiss critical chemical potential 
fic/T = 7r/3. For T > 1.5 T c , nq/T 3 as a function of [li/T exhibits an interesting 
scaling behavior with T. Right: QCD phase diagram for Nf = 4 at ima gina ry and 
real fj,B- Diamonds represent a polynomial fit to the lattice QCD data in [Z| in the 
sector of imaginary chemical potential and its analytic continuation, while the solid 
curve is the QPM result for the parametrization of the pseudo-critical line from [13| . 
See text for details of the phase structure in the imaginary sector. 



IV. SUMMARY 

In summary we show that the quasiparticle model with the chosen dispersion relations for quarks 
and gluons accounts fairly well for the quark mass dependence in the QCD equation of state near T c 
as delivered by selected lattice QCD results. We emphasize further that lattice QCD calculations at 
imaginary chemical potential, avoiding the sign problem of the fermion determinant, give valuable 
information which can be continued to real chemical potential within our model. The feasibility 
of exploring larger net baryon densities is relevant for future heavy-ion experiments at FAIR. 
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